Pair-breaking effect on mesoscopic persistent currents 
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We consider the contribution of superconducting fluctuations to the mesoscopic persistent current 
(PC) of an ensemble of normal metallic rings, made of a superconducting material whose low bare 
transition temperature T" is much smaller than the Thouless energy Ec- The effect of pair breaking 
is introduced via the example of magnetic impurities. We find that over a rather broad range of 
pair-breaking strength h/rs, such that < H/ts < Ec, the superconducting transition temperature 
is normalized down to minute values or zero while the PC is hardly affected. This may provide an 
explanation for the magnitude of the average PC's in copper and gold, as well as a way to determine 
their T^'s. The dependence of the current and the dominant superconducting fluctuations on EcTs 
and on the ratio between Ec and the temperature is analyzed. The measured PC's in copper (gold) 
correspond to Tc of a few (a fraction of) mK. 

PACS numbers: 74.78.Na, 73.23.Ra, 74.40.+k, 74.25.Ha 



I. INTRODUCTION 

Equilibrium persistent currents (PC's), flowing in nor- 
mal mesoscopic metallic rings, have been a challenge for 
both experimentalists and theorists. The persistent cur- 
rent is a manifestation of the Aharonov-Bohm effect: it 
appears when the ring is threaded by a magnetic flux, 
and it is periodic in the flux enclosed in the ringiii^ Due 
to energy-averaging and phase-coherence hmitations, one 
expects to monitor in experiment only the lowest har- 
monics in the flux quantum h/e. 

Surprisingly enough, the magnitudes of the PC's mea- 
sured on huge collections of rings (10'' copper rings'^ and 
10^ silver rings^) turned out to be larger than those ex- 
pected theoretically. The periodicity observed in these 
large ensembles is h/2e, i.e., half of the magnetic flux 
quantum. On the other hand, measurements on a single 
ring^ or on a small number— of gold rings showed the 
h/e periodicity. In the collection of 30 gold ringsi both 
the the h/2e harmonic and the h/e one were observed. 
Overall, the sign of the amplitude of the h/2e harmonic 
measured on metallic rings seems to indicate that the 
low-flux response is diamagnetid^i^ 

In the experiments on ensembles of rings,— i^i^ the av- 
erage PC was found by measuring the magnetic moment 
produced by all rings, which was then divided by the 
number of rings, N, to yield the net average current of a 
single ring. In most of the expcriments^i^i^ii the magni- 
tude of the average PC, at low temperatures, is roughly 
of the order of eEc/h. Here Ec = HD/L^ is the Thouless 
energy, L is the circumference of the ring and D — w^Zoi/S 
is the diffusion coefficient, where ^ci is the elastic mean 
free path, and vp '^s the Fermi velocity. (We consider the 
diffusive, L » Zei, case.) 

The first theoretical studies of the PC have been car- 
ried out on grand-canonical systems of non-interacting 
electronsj^i^ In these theories, the current in each ring 



is h/e periodic. The sign and magnitude of the PC of 
the individual rings vary randomly due to their high sen- 
sitivity to the disorder and to the system's size. This 
results in a very small average PC, which is dominated 
by the exponential factor exp(— L/24i)- Hence, the typ- 
ical magnitude of the current is predicted to be VN 
times the standard deviation of the PC of non- interacting 
electrons, which at low temperatures is of the order of 
eEc/h. Consequently, the persistent current carried by 
non-interacting electrons is too small to explain the large- 
ensemble experiments. Similarly, the PC predicted for 
non-interacting electrons in the canonical ensemble^ is 
substantially too small to explain the observed amplitude 
of the h/2e harmonic. 



The theory for interacting electron system: 
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pre- 
dicts h/2e periodicity of the interaction-dependent part 
of the PC. According to this theory, the average mag- 
nitude of the PC per ring due to interactions is inde- 
pendent of the number of rings. The total measured PC, 
divided by N , is thus expected to have an TV-independent 
contribution due to interactions, and an interaction- 
independent contribution proportional to iV~^/^. The 
presence of the h/e harmonic in the measurements per- 
formed on a single ring^i^ and on a fewi rings, and its ab- 
sence in large ensembles;^ are in agreement with these 
theoretical predictions. Experiments on a single ringi^ 
and on a large ensemblei^ of semiconducting rings show 
the h/e and the h/2e periodicities respectively, consistent 
with the arguments given above. 

Notwithstanding the order of the harmonics, their am- 
plitudes, in particular that of the h/2e one, remained un- 
explained for the large-ensemble measurements. On the 
other hand, the magnitudes of the h/e harmonic mea- 
sured in Refs. and|3 agree roughly with the prediction 
for non-interacting electrons, while the PC measured by 
Chandrasekhar et al^ turns out to be much larger, how- 
ever. 



2 



Here we study the PC of large ensembles, focusing on 
the role of electronic interactions. These, attractive and 
repulsive of reasonable strengths, give rise to compara- 
ble magnitudes of the averaged PC (within an order of 
magnitude), but predict opposite signs. Whereas repul- 
sive electron-electron interactionsii result in a paramag- 
netic response at small magnetic fluxes, attractive inter- 
actions yield a diamagnetic responsCfi^ as indeed seems 
to be indicated in the experiments. The magnitude of 
the PC predicted for electrons which interact repulsively 
is smaller— by a factor of about five than, e.g., the mag- 
nitude of the PC measured in copper^ The effective 
coupling strength of repulsive interactions decreases as 
the temperature decreases, due to interactions mediated 
by states whose energies are large compared with the 
temperature ] ^^1^^ This "downwards" renormalization is 
the reason for the disagreement between the theory for 
electrons interacting repulsively and the experiments..^^ 
On the other hand, the attractive interaction is normal- 
ized "upwards" at low temperatures and eventually 
leads to a superconducting state. One expects the mag- 
nitude of the averaged PC due to attractive interactions, 
i.e., due to superconducting fluctuations^ to increase 
with the strength of the interaction, or alternatively, to 
decrease as the (superconducting) transition temperature 
is reduced. Since the transition temperatures of metals 
such as copper, gold, and silver - on which the PC has 
been measured - are expected^^ to be extremely small or 
zero, Ambegaokar and Eckeri>i2i have employed in their 
estimates small values of the attractive coupling. Con- 
sequently, they came up with a magnitude for the PC 
which is again smaller by a factor of order five than the 
measured one*^ 

In order to reconcile the relatively large interaction 
required to fit the experiments with the apparent ab- 
sence of a superconducting transition, we propose that 
the rings (of e.g., copper) contain a tiny amount of mag- 
netic impurities. We show that a small concentration of 
these pair-breakers may suffice to hinder the appearance 
of superconductivity, while hardly affecting the magni- 
tude of the PC. Indeed, it seems that a small amount 
of magnetic impurities is almost unavoidable in metals 
like copper. This is suggested by recent experiments)^ 
aimed to measure the temperature dependence of the dc- 
phasing time in noble metal samples. Theoretically, one 
expects^ii^S this rate to vanish as the temperature goes 
to zero. However, it was found that the dephasing time 
may cease to increase below a certain temperature. This 
finding was attributedi^ to the presence of a small con- 
centration of magnetic impurities, which was reported to 
exist in these samples. 

As is well-known, magnetic impurities act as pair 
breakers, leading to the vanishing of the transition tem- 
perature Tc once the spin-scattering rate I/ts is larger 
than the hare transition temperature of the material with- 
out the magnetic impurities, At the same time, 
superconducting fluctuations can result in a significant 
PC provided that the lifetime of a Cooper-pair (^ Ts 



at low temperatures) is larger than the time it takes it 
to encircle the ring, ~ hjE^. (In the experiments^ii 
Ec ^ 10 mK.) Therefore, the observation that the PC is 
almost unaffected by magnetic impurities while van- 
ishes holds in the range 

n < < , (1) 

(from now on we use units in which h = 1). 

It is instructive to write the above condition in terms 
of lengths, for which Eq. ([T]) reads 

L<Ls< m , (2) 

where 

- (i?r,)i/2, and C(0) = {D/T^Y^' . (3) 

Here the magnetic-impurities scattering length Ls is the 
distance a diffusing electron covers during the time in- 
terval Ts- The bulk superconducting coherence length, in 
the absence of magnetic impurities ^(0), is the character- 
istic distance between two electrons forming a Cooper- 
pair. At low temperatures a Cooper-pair fluctuation can 
propagate a distance of the order of Lg until it is de- 
stroyed due to the scattering by magnetic impurities. 
When L < Lg the pairs are sensitive to the Aharonov- 
Bohm flux and consequently contribute significantly to 
the PC. When pair breaking occurs on scales smaller than 
the characteristics distance between two paired electrons, 
i.e.; when ^(0) > L^, then the bulk material would not 
become a superconductor. Therefore, rings made of al- 
loys which are not superconducting in the bulk due to 
pair breakers, will have PC's due to Cooper-pair fluctu- 
ations provided that Eq. ([2]) is satisfied. We show that 
the measured amplitude of the h/2e harmonic in copper— 
and gold^ rings can be understood theoretically, assum- 
ing a minute, less than one part per million, concentra- 
tion of pair breakers. Similar amounts of magnetic im- 
purities were obtained for the most purified copper and 
gold samples in Ref. [l9|. We point out that according to 
our considerations, the measurement of the PC provides 
a way to estimate T°, which may well be unreachable by 
direct experiments. 

This paper is organized as follows. In Sec. [H] and Ap- 
pendix 1^ we derive the expression for the PC due to 
superconducting fluctuations, taking into account the ef- 
fect of pair breakers. In Sec. lIIII we characterize the dom- 
inant Matsubara frequencies and wave numbers that con- 
tribute to the PC, and discuss the significant harmonics. 
In Sec. IIVI we expand the expression for the PC in the 
limits of high and low temperatures. The effect of pair 
breaking on the renormalization of the attractive inter- 
action is discussed in Sec. |Vl In Sec. IVII we present a 
detailed comparison of our results with the experimental 
data, and estimate for copper and gold. Finally, the 
results are summarized in Sec. IVIII 

In our analysis the effect of pair-breaking is brought 
about by the presence of magnetic impurities, disregard- 
ing the Kondo screening of the spins. Obviously one may 
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consider other pair breakers, such as two-level systemsr^ 
inelastic scattering;^ or magnetic fieldsi^ Other effects 
of magnetic impurities have previously been considered 
in Ref . M- 

It was suggested by Kravtsov and Altshuler— that the 
measured currents have a different source than the equi- 
librium PC discussed so far. A non-equilibrium noise, for 
example, a stray ac electric field, can cause a dc current 
by a rectification effect. In Ref. [2y it was shown that the 
measured signal^ may be explained provided that there 
exists such a non-equilibrium noise. This mechanism is 
different than the one suggested by us. 



II. DERIVATION OF THE PERSISTENT 
CURRENT 

The PC is obtained by differentiating the free energy 
of electrons residing in a ring with respect to the mag- 
netic flux enclosed in that ring. In this section we derive 
the term in the free energy which results from super- 
conducting fluctuations. The system consists of diffusing 
electrons which interact with each other attractively, and 
are also scattered by magnetic impurities that couple to 
their spin degrees of freedom. We use the Hamiltonian^i 



H 



drUt(r) (7^o + "i(r))(5„^+U2(r)S.^°T ^Jr) 
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(4) 



in which the last term represents the attractive interac- 
tion, with coupling g (> 0). The spin components are a 
and 7, and <t is the vector of the Pauli matrices. The 
free, spin-independent, part of the Hamiltonian is 



Ha {-iV - [27r/L]0x)V2m - ^ 



(5) 



where m is the electron mass, fi is the chemical potential, 
and (j) is the magnetic flux through the ring, in units of 
h/e. The unit vector x points along the circumference 



of the ring in the anti-clockwise direction. The scatter- 
ing, by both nonmagnetic and magnetic ions, is assumed 
to result from Ni point-like randomly-located impurities, 
such that 

iti(r) + U2(r)S • cr 

Ni 

^ {S{r -^^)-y) (^1 + "2Sr, • <t) , (6) 

i=l 

where V is the system volume. 

We calculate the partition function Z using the method 
of Feynman path integrals combined with the Grass- 
mann algebra of many-body fermionic coherent states^ 
in which the superconducting order-parameter is intro- 
duced by the Hubbard-Stratonovich transformation!^ 
Details of this procedure are given in Appendix [Xj As is 
shown there, the partition function is (the temperature 
is denoted by T) 



z = Zon(i-f^n(q,^) 



where the polarization^ 

n(q, ^ X! '^ajKcjajifi, v) , 

consists of the Cooperon-dominated contributions 

K^a-i{(\, ^) = 'Y {Gaa' (ki + q, k2 + q, + z/) 
ki,k2 

X eQ'^'G^^'(-ki, -k2, -t^)) . 



(7) 



(8) 



(9) 



Here e is the anti-symmetric tensor, Saa — 0, and e-\i ~ 
— Ejl = 1, and G denotes the particle Green function. 

In Ref. [U the polarization n(q = {),v = Q) was calcu- 
lated from the Dyson equation for the Cooperon. Their 
calculation can be extended to general q, v 



= Gc,a(k + q, w + i/)G^^(-k, -u) [eo,^ + iYi(wi(5„a' + "28 • o-««')(wi^77' + ^*2S • cr77')i^„^,^, (q, v)] . (10) 

k 



Here, Gq^ = Sa^Gaa is the Green function averaged over 
the impurity disorder and spin components (which makes 
it diagonal in spin-space). Averaging over the impurity 
spins. 



(Ul^Q 



-U2S • <T««')(ui 5^^' + U2S • 0-77') 



(11) 



is carried out employing Si 
l)/3 (where «, j = x,y,z). 



and SiSj ~ Sij S{S 



Following Ref. [2l| we assume that 



^cfy-^uj^ and 
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then using cr""'o'J''' ^a'-,' = - 
n(q,i.) =^/C(q,i^) , 



-if-a-i wc obtain 



K^(ci,v) = [1 + (27rAA(0)r_)-iA'^(q,z/)] 

X ^G(k + q,cj + i^)G(-k,-cj) , (12) 

k 

where the averaged Green function is 

G(p, w) = [icj - (pVSto - + isgn(tj)/2r+]"i . (13) 

(The spin indices are suppressed since G is independent 
of them.) In Eqs. ^ and (US]), 

— = 2T^M{^)Ni{u\ ± S'(S' + \)u\) , (14) 

where A/'(0) is the extensive density of states at the Fermi 
level. (Note that t+ is the elastic mean free time.) Using 
Eq. (fT5|) to calculate the sum over k in Eq. yields 

^ G(k' + q, + i^)G(-k', -a>) = 27rA/'(0)r+ 



X 6I[w(lj + i/)](l - T+|2tj + j/| - £'q^r+) 



(15) 



Upon inserting this expression into Eq. (jl2p and solving 
it one finds 



K^{ci, v) = 27rAA(0)6'[cj(tj + z^)] 
X (i?q2 + |2Lj + z/| + 2/r«)-i , 

where 1 /r^ is the pair- breaking rate 

— = 27rAr(0)iV,S'(5' + \)u\ . 



(16) 



(17) 



When T+ ~ T_ most of the disorder is due to the non- 
magnetic part. This, together with the assumption^ 
{|2(jj -I- I'l, £)q^} ^ l/r+ were used in obtaining Eq. (fTH]) . 

The summation in Eq. over the Matsubara fre- 
quencies can be written explicitly as 



T 



Af{Q) 



n(q,^^) = £ 



71=0 '- 



"+2 + i^T 



(18) 



Note that Eq. ([T8|) includes also the negative Matsubara 
frequencies. This sum does not converge and therefore a 
cutoff is required. The cutoff frequency on the attractive 
interaction is the Debye frequency coo, and consequently 
the sum is terminated at n = luu/2t:T. As a result, the 
polarization is given by 



T 



Af{0) 



2^ 



\u\+2/Ts + Dq^ 



,^/l \v\ + 2/Ts + De 
\2 AttT 



(19) 



where is the digamma function. 

We next express the polarization in terms of the bare 
transition temperature of the system. This is the temper- 
ature at which Z/Zq diverges for = and the smallest 
possible |q|, in the absence of the pair breakers and the 
magnetic flux, 



5AA(0) b 



27rT" 



- * 



(20) 



Since ujjj 3> {T^, T} we may use the asymptotic expan- 
sion of the digamma function. 



^{x > 1) ~ ln(a;) , 



(21) 



to obtain 



Z = Z, 



n 



In 



TO J 



2/ts + Dq^ 
AttT 



-*'2 



(22) 



The effect of the pair breakers is represented by the term 
2/ts in the argument of the digamma functions. 

As is mentioned above, the persistent current is given 
by 



/ = (e/27r) dT\iiZ/d4> 



(23) 



The flux enters the expression for Z through the lon- 
gitudinal components of the momenta, see Eq. (|A14|1 . 
In our ring geometry, only momenta of zero transverse 
components contribute significantly to the current, since 
the contribution of momenta of higher transverse compo- 
nents can be shown to decay exponentially, as a function 
of the ratio of L and the transverse dimension (e.g., the 
height) of the ring. 

As is seen in Eqs. ([22]) and (|23p . the PC consists of 
two parts. The first arises from differentiating Zq and 
is the ensemble averaged PC of non-interacting, grand- 
canonical, normal metal ringsi^ This contribution is much 
too small to account for the measured amplitude of the 
h/2e harmonic [see Sec.|I], and therefore will be omitted 
in the following. The other part of the PC comes from 
the free energy due to the superconducting fluctuations. 



-2eEr 



(n + 20)^^(F(n,i/)) 
ln(r/r0) + *(nn,z.))-vl/(i) ' 



(24) 



where we have introduced the function 

1 , \v\+2/ts , ^E, 
2 



^K^) = ^ + ^4z^ + ^(" + 20)' ■ (25) 



A-kT 



T 



In particular, one notes the h/2e periodicity in the flux. 
Indeed, upon employing the Poisson summation formula 



sin(47rm0) 



„ Jo 



xsm{2'Kx)'9'{F{x, ly)) 



(26) 
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where 



F(.T,i.) = i + 



Clearly, the fluctuation-induced PC decreases as the pair- 
breaking strength increases. Our central result is that 
this decrease may be far less than the one caused in the 
transition temperature. 

In order to compare the dependence of the PC and of 
the transition temperature on the pair-breaking strength, 
wc use the expression^ for the transition temperature in 
the presence of both pair breakers and magnetic flux 



In 



TP 



1 



2ttTcTs 







(28) 



Here (p is in the range —1/2,1/2, modulo unity;^ We plot 
in Fig. [T]the amplitude of the h/2e harmonic of the PC, 
as well as the transition temperature (in the absence of 
the flux) as functions of the pair-breaking strength, using 
the dimensionless parameter 



l/7rT°7 



(29) 



The transition temperature is reduced due to pair break- 
ing and vanishes at s = 1/2je, where 7b is the Euler 
constant. In contrast, for Ec ^ 1/ts the PC is hardly 
affected for these values of pair-breaking strengths. 
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FIG. 2; The amplitude of the h/2e harmonic in units of J* = 
—eEc, as a function of the temperature, for two values of 
T^/Ec and several values of s. Note that the s = curve in 
the upper panel is valid only for T/Tc > 1 + Gi, where Gi is 
the Ginzburg parameter. 



III. THE DOMINANT FLUCTUATIONS 

Our result for the PC [see Eq. (|24)) ] consists of inflnite 
sums over the frequencies and over the momenta. One 
naturally asks oneself whether the characteristic features 
of the expression are not given by the first few mem- 
bers of each sum, notably the static, = 0, regime. 
It turns out that this is not the case over most of the 
relevant range: to obtain the correct magnitude of the 
fluctuation-induced PC, numerous frequencies and mo- 
menta are required. 

In order to study this aspect, it is convenient to express 
the PC in a form which is more amenable to numerical 
computations. To this end we write Eq. (|26p as 



FIG. 1: The h/2e harmonic (full line) and Tc/T° (dashed 
line) as functions of the pair-breaking strength, displayed on 
a logarithmic scale. The current, in units of I(s = 0), is 
plotted for T = Ec and T° = O.lEc. The PC reduction at 
s = 10 corresponds to 1/ts = ttEc. 

Figure [2] portrays the PC plotted by numerically eval- 
uating Eq. (|26p . In each of the panels the upper curve is 
drawn for s — 0, while the second curve corresponds to 
a pair-breaking strength [see Eq. (pS)) and Fig. [T] which 
is large enough to destroy Tc. Nonetheless, the PC is 
hardly affected as long as > L [see Eqs. ([2]) and 
The considerably-reduced PC due to a small Lg is pre- 
sented by the dash-dotted curves, which correspond to 
Lg ~ 0.5L. The effect of the temperature on the mag- 
nitude of the PC is manifested by its dependence of the 
ratio L/Lt, where Lt is the thermal length. 



Lt = VD/T , 
or equivalently the ratio T/Ec, sec Fig. [51 



(30) 



^ 2ieT ^ . r , 

I = > sm(47rm0) >^ / axe 

m=l u 



IT 



X — In [^{F{x, v)) - HT^JA^eT)] , (31) 

where the function F is given in Eq. (|27p . The x- 
integration is carried out by closing the integral in the 
upper half of the complex plane. Two sets of simple poles 
can be identifled in the integrand of Eq. ((3T|) . These sets 
result from (a) the zeros and (b) the poles of the argu- 
ment of the logarithmi^i The first set of poles, denoted 
by a:fe„, is given by 



*(FiJ = in(r,747^,r) . 



(32) 



The second set consists of the poles of the digamma func- 
tion. These are denoted by x^^^y^, and are obtained from 
the relation 



F' 



£=0,1,2, 



(33) 
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The index ^ runs over the poles in each set. The two sets 
•-"^ ^poic/zoro gi'^*5n by Eqs. ((5^ and ([55]) . are shown in 

Fig.i' 
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FIG. 3: The digamma function (solid line) and ln(ri'/47BT) 
for T° /T = 0.6 (dashed line). The first three solutions -F^ero of 
Eq. p2|l are marked on the x axis with their indices indicated 
below it. The first values of the set Fp^i^, Eq. (|33p . are marked 
by arrows. 



Performing the Cauchy integration, the current takes 
the form2^ 



/ = - 4eT ^ sin(47rTO0) 

m— 1 
oo 

X ^ ^ [exp(27ria::f^„) - exp(27rza;^„ij] . (34) 



Here x^^^^^^^^.^ depends on the Matsubara frequency and 
the harmonic index m, 



polc/zcro 



T 



2t:E, 



1 



2ttT 



2R 



polc/z 



1/2 



(35) 



Note that all the exponents (27ria;) in the two series in 
Eq. ([34)) are negative, and their absolute value increases 
with increasing v, I or m. As can be seen from Fig. [31 
for each pair of poles F^ro > ^poiei consequently 
I^^Liol < l^^poicl- This ensures that the term in the square 
brackets of Eq. ([34)) is positive, and hence the response 
of the ring to a small flux is diamagnetic, as it should be. 



A. The dominant imaginary time fluctuations 



The dominant terms in Eq. ([34]) are those for which 
the absolute value of x is smaller than unity, but if the 
absolute values of all x are larger than one only the 
smallest [xf^o(i^ = 0, m — 1)] is the dominant one. 
The absolute value of the exponents (which are given by 
27r|a;zei.o/poiel) is at least (|z^|/i?c)^^^. Thus, the frequen- 
cies that contribute mostly to the current are those for 
which < lOEc- The proportionality factor, of order 
10, had been determined numerically and resulted from 
the square-root structure of the exponents, see Eq. ([55)) . 
At high temperatures T > Ec the system is dominated 



by the classical fluctuations - namely, by the first (low- 
est energy), 1^ = 0, Matsubara frequency. The effect of 
the quantum fluctuations for which i/ ^ increases as 
the temperature decreases. This tendency has an excep- 
tion in two cases. First, for very strong pair breaking 
1/ts > {T,Ec,T^/Ec} the significant quantum fiuctua- 
tions that have a dominant contribution to the PC are 
bounded by < \/~E^tI. Second, in the case of small 
or zero pair breaking when T ^ Tc, only v = Q la the 
dominant frequencyi^ 

When Tc is finite, the n = v = pole of the partition 
function, Eq. ([7]), is the most dominant one as T ^ Tc. 
Consequently, in this low-temperature regime physical 
properties, including the PC, are determined only by the 
v = Q fluctuations, pertaining to the static Ginzburg- 
Landau free energy. We find however, that in the case of a 
vanishing Tc, quantum fluctuations, for which v ^ {] have 
a significant contribution to the PC at low temperatures. 
Indeed, the quantum fluctuations of a system with no 
magnetic impurities and for which > /{IGttjeEc) 
have been recently invoked in the context of the "strong" 
Little-Parks oscillations, see Ref. [H- 



B. The dominant spatial fluctuations 

High Matsubara frequencies involve many spatial fre- 
quencies q. Thus, at low temperatures and for a van- 
ishing Tc, many wave vectors contribute to the PC. We 
have estimated numerically their number by comparing 
the PC computed with a relatively small number of fre- 
quencies and wave vectors with the exact result Eq. ([34)) 
for T = T° = 0.1£^c and s = 1. In this case we have 
found that ^ 100 Matsubara frequencies are required. 
The highest momenta, Eq. (|A14[) . that contribute sig- 
nificantly are given by |n| ^ (1, 5, 100, 1000) for the fre- 
quencies i'/{2itT) = (0,5,10,100), respectively. Figure 
ID shows the PC as computed from Eq. ([M)) for different 
maximal |q| values and without limiting the range of v. 
It is thus seen that in the whole range of 4> the persistent 
current is not mainly determined by the lowest momenta, 
even when the size of the system L is smaller than the 
thermal length Lt, Eq. ([30)) . This is different than the 
situation in the calculations of other properties (for ex- 
ample, weak-localization corrections^ to the conductiv- 
ity), in which L <^ Lt is taken as a sufficient condition 
for using only q — 0. We point out however that the 
slope of the PC at = appears to be describable using 
the smallest wave number only. 



C. The dominant harmonics 

Examining the series in Eq. ([M)) one can see that 
the maximal harmonic of the fiux, mmax, that has still 
a significant contribution to the current is given by 
min{ Ec/T, ^/E^} or by one if the first two values 
are smaller than unity. This condition can be expressed 
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l/eE^ L - -j:, 

c 



-0.2 



-0.1 







0.1 



0.2 



FIG. 4: The PC as computed from Eq. ([241, with the sum- 
mation over n cut at 1000, 3, 1 (solid, dashed and dash-dotted 
curves, respectively). The plots are for T = = O.lEc and 
s = 1. 



in terms of lengths by 



minjL^/i, Lf/L} , or 1. 



(36) 



The upper limit on the harmonics results from the fact 
that the m'th harmonic is associated with paths that 
encircle the ring (coherently) m times and hence their 
length is at least mL.— The sinusoidal shape / (x 
sin(47r(/)) at high temperatures is modified due to higher 
harmonics as the temperature decreases. In the absence 
of magnetic impurities (upper panel in Fig. the low- 
temperature current as a function of the flux attains a 
sawtooth shape. Such a behavior is predicted also for the 
equilibrium PC in superconductors at zero temperature*^ 
and for the persistent current in a clean system of non- 
interacting electronsi^ In the presence of pair-breakers 
the upper bound on the harmonics Eq. ()36|) prevents the 
current from reaching the sharp sawtooth shape. This 
suggests, in principle, a way to experimentally confirm 
the role of pair breaking for this problem. In the lower 
panel of Fig. [S] the current of a system with L ~ is 
plotted for several temperatures. At temperatures below 
O.lEc the shape of the current does not change anymore. 



IV. THE TEMPERATURE DEPENDENCE 

Here we study the PC in the low and high-temperature 
regimes. In particular we find that the PC decays expo- 
nentially as the length of the ring exceeds the thermal 
length Lt or the magnetic impurity scattering length Ls , 
whichever is shorter. 



A. High-temperature regime, T max{l/rs, T^, iSc} 

When the temperature is much higher than all rele- 
vant energy scales, i.e., T 3> max{l/rs, T^, Ec}, the lead- 
ing contribution to the double sum in Eq. ([M]) comes 
solely from the first pole x^^^o of t^^^ lowest Matsubara 
frequency, i/ = [see Eq. ([55]) ]. In this temperature range 
the h/2e harmonic, corresponding to m = 1, is the dom- 
inant one. 

As the temperature increases, the horizontal line in 
Fig. [3] representing ln(T^/47£;r) moves further down, so 



l/eE 



l/eE 



1/x =0 

s 



\ 



*' * - - 



1/7tT°t =1 



-0.2 



-0.1 







0.1 



0.2 



FIG. 5: The current, in units of eE^, as a function of the flux 
for /Ec — 0.1 and for several temperatures, T/Ec = 5, 1 
and 0.15 in the solid, dashed and dash-dotted curves respec- 
tively. In the lower panel the dash-dotted curve corresponds 
to T/Ec = 0.1. For s = the current attains the sawtooth 
form (upper panel) which is lost for s = 1 (lower panel). 



that approaches zero. We use the expansion of the 

digamma function for small arguments in Eq. (j32p and 
obtain 



[in(r,"/47|r) 



(37) 



Upon substituting this result in the dominant term of 
Eq. ([M]) . we obtain the current in the form 



4eTsin(47r(?!)) 
L 
Lt 



I - 



X expl — -— 



2tt 



7rln(47|T/rO) 



1/2 



(38) 



We compare the full result, Eq. with the high- 

temperature approximation Eq. psp in Fig. [B) The dif- 
ference between the contributions of the first xf^ol^ — 0) 
and the second x^i^{i' = 0) poles to the PC is the ab- 
sence of the third term, which includes a logarithm [see 
Eq. ([55]) ]. in the exponent of the latter. Therefore, this 
approximation improves as increases. 




FIG. 6: The amplitude of the h/2e harmonic is plotted in 
units of I* — —eEc as a function of the temperature, for 
/ Ec — 0.1 and s = 1. The exact results can be approxi- 
mated by Eq. (f38)l for T > Er. 
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B. Low-temperature regime, Tc <^ T <S {1/ts,Ec} 

In the low temperature regime the argument {F) of 
the digamma function and its derivative is much larger 
than unity [see Eq. ([?/)) ] , so that we can use their asymp- 
totic expansions ln(F) and respectively. Substitut- 
ing these approximations in Eq. (1261) gives 



E 



iT sin(47rm(^) 

X sm{2'Kx)dx 



m—l 



In 



(39) 



where a^,^ = m^{\v\ + + 2iTT)/{An'^Ec). For T > 
Tc the denominator in Eq. ((5^ docs not vanish. Then 
the term x'^ + am,^ in the logarithm in Eq. ((39)) can be 
replaced by aam,i^, with, say, 1 < a < 3. Consequently, 



/In 



4ery~^ sin(47rm0) e ' 



TP 



27rT 



(40) 



Since T ^ Ec the summation over v can be replaced by 
an integration. Approximating again the logarithm by 
its value at the dominant ly of the integration, yields 



77 ^ ^ 



sin(47rm(? 



/27ri2 

Lrp 



2L2 



X e 



/In 



(41) 



where ^ = max{l, 2m2/rs£'c}. 

We compare in Fig. [7] the low-tcmpcrature approxi- 
mation, Eq. dill), with the full result, Eq. ([34]). As one 
can see from this comparison, the flux dependence of the 
PC as well as its amplitude are well approximated by 
Eq. gl]). 



l/eE 



0.5 

" 
-0.5 





exact 
appr. 









-0.2 



-0.1 







0.1 



0.2 



FIG. 7: The current in units of eEa as a function of the 
magnetic flux 4>, plotted for T = CIT^ = O.OlSc and s = 
1. The low-temperature approximation Eq. (|4ip is compared 
with the exact result Eq. p4|) . We take a = 3 in the logarithm 
of Eq. dHJ. 



V. RENORMALIZATION OF THE EFFECTIVE 
INTERACTION 

In this section we calculate the PC to first order in 
the interaction, in order to see whether it suffices to ex- 
plain our full result. To first order in the interaction, the 
contribution of superconducting fluctuations to the free 
energy [see Eq. ^\ is 



(42) 



The PC resulting from Eq. ([IS]) has the same form as 
Eq. ((24|) , except that the denominator in the latter is re- 
placed by the bare interaction gA/'(0)/y. Had we tried 
to to fit the experimental data of Refs. and us- 
ing Eq. (|42)) . we should have taken the implausible ra- 
tio Ec ^ O.Iwd [see Eq. (|20p ]. This first-order approxi- 
mation fails because of screening effects, which increase 
the magnitude of the effective attractive interaction as 
the temperature decreases. Very roughly, the renormal- 
ization of a dimensionless interaction A, from a higher 
frequency scale a;> to a lower one, tt'<, is given byi^ 



A(c^<) 



A-^(^>)-ln^ 



(43) 



For attractive interactions A is positive and the high fre- 
quency scale YA lod- At T = and 1/ts = 0, the attrac- 
tive interaction should diverge. Using this to eliminate 
X{ujd) (= gJV{0)/V), we obtain that for T^" < < ivd, 



A(c.)^l/ln(^/r,") 



(44) 



Replacing in the first-order approximation for the current 
the bare interaction by the effective one, Eq. (UH), gives 



list — — 



8eE,. 



E 



sin(47rm0) 



xV/ dxxsm{2TTx)^'{F{x,iy)) . (45) 

The effective interaction is renormalized upwards with 
decreasing energy and, for the bulk and no pair breaking, 
it blows up at T°. For I/t, > r°, this renormalization 
stops at 1/ts and Tc disappears. In the mesoscopic range, 
the Thouless energy, Ec, becomes a relevant scale and it 
may be expected (as is borne out by our results) that the 
PC at low temperatures is determined by the interaction 
on that scale, as long as Ec > 1/ts. Once l/r, > Ec, 
we expect the renormalization to "stop at l/r^" and the 
PC to be depressed. Thus, the relevant range for our 
considerations is < 1/ts < Ec- Using these bounds 
on the energy scale of the renormalized interaction in the 
first order calculation Eq. (|45p . gives a good agreement 
with our result Eq. (|M|) . In Fig. [5] we plot the amplitude 
of the h/2e harmonic as a function of T/Ec, calculated 
from the full expression Eq. (j26p (thin curves) and from 
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T/E 



FIG. 8: The first-order approximation for tiie h/2e iiarmonic 
of the current Eq. (|45|l (bold Unes) is compared with the exact 
result (thin lines). Here — Q.lEa- In drawing the former, 
we have used the simplest expression for the cutoff w = T + 

+ I/ts . 



the first-order approximation Eq. (P5|) (bold curves). The 
plotted curves are for = O.lEc- 

A more precise expression for the renormalized attrac- 
tive interaction depends on q, u of the order-parameter 
fluctuation. The renormalized attractive interaction 
A(q, v), obtained from an infinite series of diagrams con- 
taining Cooperon contributions, is given bj*^ 



A(q,i.)= X~\ujd) 



T 



AA(0) 



(46) 



Upon substituting Eq. p9|) in Eq. (|46p one can identify 
A(q, ly) from our result, e.g., by comparing Eq. (j45p with 
Eq. jSl). 



VI. COMPARISON WITH EXPERIMENTS 



Theoretically, only static magnetic fields have been 
considered here. However, experiments have been car- 
ried out with an ac magnetic field. In the experiments 
on copper— and gold^ the sweeping frequencies of the 
magnetic field were very low (0.3 Hz and 2Hz respec- 
tively). Thus, one expects that the measured PC could 
be explained using a theory for a static magnetic field. 
In the experiment on silver, on the other hand, a very 
high sweeping frequency of the magnetic field was used 
(217 MHz). It is plausible that in order to explain the 
results of Ref. one may not confine oneself to a static 
magnetic field. We therefore do not attempt to explain 
the experiment of Ref. 0. 

Here we explain the h/2e signal observed in copper— 
and gold^ using our result Eq. ([34|) . In the left six 
columns of Table [T] we summarize the experimental pa- 
rameters for the (/i/2e)-periodic signal4S 







T 


I/eE, 


L 


Lip 


min 


Coppes^ 


15 mK 


7 mK 


1 


2.2 ^m 


2 fj,m{1.5 K) 


a few mK 


Goldi 


4.9 mK 


5.5 mK 


0.65 


8.0 ^im 


16 /im(0.5 K) 


a fraction of a mK 



TABLE I: Experimental parameters in the left six columns. The magnitude of the h/2e periodic current (column 4) is given for 
the lowest temperature (column 3) reached in the experiment. The dephasing length is given together with the temperature 
at which it was measured. The last column is our estimate for a lower bound on according to Eq. (|34|l . see also Fig. [9] 



The metals used in the experiments are not supercon- 
ductors at any measured temperature in their bulk form. 
Therefore, it is not possible to obtain theoretically a large 
enough PC (to match the measurments^) due to the at- 
tractive interaction without pair breaking: the required 
Tc ^ 1 mK are too high to be considered as realistic. We 
suggest that the bare transition temperature may indeed 
be of the order of a mK, but the transition tempera- 
ture of the real material is considerably reduced due to 
pair breakers. Together with this assumption, the neces- 
sary condition to fit the experiments is l/r, > '!rT^/2jE 
so that Tc vanishes or is very strongly depressed;^ see 
Eq. This condition can also be written as 



Note that we need > L in order not to depress the 
PC [Eq. ([34]) ]. The upper limit on T°, corresponding 
to the equality in Eq. (j47|) . is given by the solid line 
in Fig. [HI The values for T^/Ec that correspond to a 
vanishing Tc are in the region below this line. In the 
dashed and dash-dotted curves in Fig. [9] different values 
of L/Lg are matched with an appropriate T^/Ec so that 
the measured values in columns 2 — 4 of Tabic [J remain 
the same. 

The monotonically increasing shape of the curves in 
Fig. [9] results from the fact that higher values of r° / Ec 
are required to describe the experiments as L/Ls in- 
creases. The minimal T^'s correspond to the points 
where the dashed and the dash-dotted lines cross the 
solid line. In this way we obtain estimates of the lower 
bounds on the value of T° for copper and gold. These 
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L/L 

s 

FIG. 9: The bare transition temperatures corresponding to 
the measured PC as a function of L/Ls. The dashed and 
dash-dotted curves correspond to the PC measured in copper 
and gold, respectively. The solid curve gives the maximal 
possible satisfying = 0. 

lower bounds arc given in the seventh column of Table [H 
These estimates of T° are very sensitive [see Eqs. (PT|) 
and ([45|) ] to the experimental parameters. For example, 
Ref. points out that the measured values of the PC are 
correct up to a factor of two. The exact minimal value of 
T° that satisfies Eq. (gT]) for copper, based on the values 
quoted in TableU is 4mK. However, assuming half of the 
value reported in Ref. for the PC, results in a minimal 
of about 0.3mK. The curves in Fig. [9]ignore the error 
bars in the experiments. Thus, the values of T^/E^ in 
this figure should be considered only as rough estimates. 

Besides spin-flip scattering from magnetic impurities, 
decoherence of the electrons is also caused by other pro- 
cesses, e.g., electron-phonon inelastic interactions. Hence 
Lg is always larger or of the order of the dephasing length. 
A lower bound on Lg is given by equating it to the mea- 
sured L^. Those values (see Table |I| are small enough 
to fulfil the condition in Eq. (|T7)) . In other words, we 
could account for the data of copper and gold since the 
measured L^p was small enough. This is not the case for 
silver^ where L/Lip{~ 0.3) is too small to explain the re- 
sult I{T = 4.6^;^) = l.QeEc using Eq. Our theory 
is not applicable to that experiment. We believe that the 
reason for that, as explained above, is the high frequency 
used in that experiment. 



VII. DISCUSSION 

In our result for the PC, Eq. ([M)) . there appears the 
bare transition temperature and not the one reduced by 
the pair-breaking mechanism. Therefore we propose the 
scenario, in which the bulk transition temperature van- 
ishes due to the pair-breaking mechanism, while the PC 
is dominated by a relatively high attractive interaction. 

The bulk Tc vanishes due to pair breaking for Ls < 
^(0). However, we find that the PC may still be hardly 
affected by pair breaking. The physical reason for that 
is that as long as Ls > -^^ the Cooper pair fluctuations 
can complete a circle around the ring before being mag- 
netically scattered, and hence respond to the Aharonov- 
Bohm flux. The PC is immune to pair breaking in the 



regime given by Eq. ^ where the bulk form is normal. 
This is demonstrated in Figs. [T] and [H 

In the pair-breaking regime given by Eq. ^ , the upper 
bound on the dominant quantum fluctuations {v ^ 0) 
is determined by the Thouless energy. Dominant fluc- 
tuations of high Matsubara frequencies necessitate high 
wave numbers. Therefore, at low temperature T <ti Ec 
high wave numbers are involved too in the dominant 
fluctuations [see Fig. 0] , in contrast to the effective di- 
mensional reduction occurring in other phenomena when 
L ^ Lt, notably weak-localization correctionsi^ The 
maximal number of flux harmonics that contribute to 
the PC, Eq. ([55)1 . is bounded due to thermal fluctuations 
and due to spin-flip scattering. Consequently, in a sys- 
tem with magnetic impurities, even at zero temperature 
the PC may not have the sawtooth shape, which appears 
for the PC's without pair breaking, see Fig. (5] 

The effective interaction is renormalized upwards 
with decreasing energy; For the bulk it stops at ~ 
max{T^,l/Ts} (which explains why Tc disappears for 
l/'''s ^ T'c). In the mesoscopic range, Ec > T^, the 
Thouless energy sets another bound for the energy scale 
at which the renormalization stops. In Sec.|V]it is shown 
that these considerations agree with our result for the 
PC Eq. (1211), see Fig. [H] 

We found that in the high-temperature regime, the 
PC decreases exponentially with L/Lg or with L/Lt, 
whichever is larger. The explicit exponential decay of 
the PC with L/Lg in both the high and the low temper- 
ature regimes [Eqs. (pS)) and (|¥T|) respectively] for L > Lg 
is in agreement with the qualitative argument of Eq. 
Note that Eq. (HT1) is applicable only at very low tem- 
peratures, such that T <C {T^,Ec}. The experiments 
on copper— and gold^ rings correspond to ^ ImK, 
thus Eq. (|4ip can be used only at very low temperatures 
T <C ImK. In the experiments the lowest temperature 
was ~ 10 mK, and therefore the measured PC cannot be 
precisely fltted by the approximate expression Eq. (|4ip . 
In the low-temperature regime the dependence of the PC 
on is logarithmically weak [see Eq. (|4T|) ]. This weak 
dependence explains why in Ref. [Tol . where the transition 
temperature was taken as lO^K (in the absence of pair 
breaking), the result was smaller only by a factor of ~ 5 
compared with the experiment 

Interestingly enough, it follows from our work that by 
measuring the PC and the pair-breaking strength one 
may determine , which would be directly measurable 
only if enough low-temperature pair breaking could be 
eliminated. This elimination is very hard to achieve in 
some materials. Our result Eq. (|26p can explain the large 
PC of Refs. i and 0, with Ls value larger than (or of 
the order of) the measured L^p (see Table |T] and Fig. [S]). 
Even though Lg was not measured in the PC experi- 
ments, we obtain a lower bound on the bare transition 
temperatures for copper and gold. These minimal T^'s 
correspond to minimal pair-breaking strength given by 
Ls ^ 5 /xm in the copper sample^ and Lg ^ 25 ^m in 
the gold samplcfi The fltted maximal Ls's can be caused 
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by a very low (less than one part per million) concentra- 
tion of magnetic impurities. These concentrations seem 
appropriate for the purest copper and gold samples avail- 
able experimentallyfiS Although, a full consideration of 
the effect of the magnetic impurities, including Kondo 
physics, is still necessary. 

Our result concerning the fundamentally different sen- 
sitivities of Tc and PC's to pair breaking is valid regard- 
less of the situation in specific materials. Our idea can be 
tested, for example, by measuring the persistent currents 
in very small rings made of a superconducting material 
whose transition temperature is known, as functions of 
possible pair-breaking mechanisms. For Ec > 100 mK, 
say, and a material with T° of a few 10 mK, the range 
of pair breaking which satisfies Eq. ([1]) becomes easier to 
control experimentally. 
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APPENDIX A: DERIVATION OF THE 
PARTITION FUNCTION 

Here we derive, using the method of Feynman path 
integral, the partition function, Eq. ((T]). In terms of 



the Grassmann variables ■0a(r, r) ['^^(r, r)], the parti- 
tion function reads 



Z= / Z?(V'(r,T),^(r,r))exp(-5) 



(Al) 



where the action S is 



J dr j\T[iPa{r,T)dr^4r,T)+nir,T)] . (A2) 

Here (3 ~ 1/T and Ti is given by the integrand of the 
Hamiltonian, Eq. (|4]), with Grassmann variables (of the 
same imaginary time) replacing the creation and annihi- 
lation operators. Introducing the bosonic fields A(r, r) 
via the Hubbard- Stratonovich transformation, the parti- 
tion function takes the form 



Z = j i?(V'(r,r),^(r,T))i?(A(r,r),A*(r,T)) 
X exp(-5) , (A3) 

where the differential of the bosonic field A(r,T), 
i:'(A(r,T), A*(r,r)), contains a factor of l3V/T:g. The 
action S is given by 



-i4'(r,T)G-J;,_,*(r,T)) , (A4) 



where ^ = (i/j|, -01, -011 "^i); and the inverse Green func- 
tion (at equal positions r and equal imaginary times 
t) is 



r— r';r— r' 



-dr-h'\ -2u2S- A 

-2u2S+ -dr-h'l -A 

-A* -dr + h~'^ 2u2S+ 

A* 2w25_ -dr + h' 



AT 



A 



Here h'^'^ = Wo(±0) + ui + sgn{a)SzU2, and S± = (Sx ± iSy)/2, where sgn(|) = 1 and sgn(j) 
The integration over the fcrmionic part in Eq. (jA3P yields 

^^,^|A(r,.)- 

We expand Trln(/3G^^) up to second order— in A 

1 



Z = J Z?(A(r,T),A*(r,r))exp(iTVln(/3G-i)-y' dr^^dr^^^^) 



Trln(/3G-i) = Trln(/3G-^) 



my 



dvdv'dTdr' Tr 

r 



Gp(r',r';r,T)A(r,r)G,(r,r;r',T')At(r',r') . (A7) 



(A5) 



(A6) 



The inverse Green function for non- interacting electrons, the right hand side of Eq. (|A7|) . which is zcroth-order in 



Gq , is given by Eq. (jASp for A = 0. The first term on 
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A, gives rise to the partition function of non- interacting 
electrons, Zg = det(/3G'(^^). 

In Eq. (jA7p . Gp (Gh) is the particle (hole) Green func- 
tion. These functions are the solutions of 

G;/,(r,r)G,/,(r,T;r',T') =<5(r-r')<5(T-T') , (A8) 

where G~J'j^ are defined in Eq. (jASp . As can be seen in 
that equation, the particle and the hole inverse Green 
functions are related to one another by 



G^i(r,T,0,5+,5_,5,) 
= -G-\r,-T,-cf>,S^,S+,S,) 



Therefore, 



(A9) 



Ghiv,T;r',T',<f>,S+,S-,S,) 
= -Gp{r,-T;r\-T\-^,S-,S+,S.) 
= -Gp{r,T';r\T,-(l),S-,S+,Sz) , 



(AlO) 



where in the last equality we have used time-translational 
invariance to shift r and r' by r + r' . Reversing the sign 
of the flux (j> together with interchanging r and r' leads 
to the relation (the superscript t denotes the transposed 
matrix) 



Gpir,T';r',r,-c^,S^,S+,S,) 
= Gl{r\T';r,T,cj,,S+,S^,S,) 



(All) 



We have used Eqs. (|A10p and (|A11|) to replace the hole 
Green function in Eq. (jA7| by a particle Green function. 
Then, in momentum representation, the second term of 
the right hand side of Eq. (|A7|) reads2£ 



Trln(/3G-i) 



2nd 



J2 J2 Tr [Gp(ki +qi,k2 + q2,t^ + i^)A(q2,i/)G[,(-ki,-k2,-tj)At(qi, 
11,12,^' ki,k2,i^ 



r 



(A12) 



The flux dependence is incorporated into the momenta p, 
where p^/2m — are the eigenvalues of Ti.o{4>). Thus, the 
longitudinal components of the momenta in the Green 
function G have the form 

2TT{n + (j))/L , (A13) 

while those of the momenta in the boson field A are 



2Tr{n + 2<P)/L , 



(A14) 



where n is an integer. The Matsubara frequencies of the 
Green functions, cu + f and — tj, are fermionic [~ 7r(2?i + 



l)r]. The order-parameter fluctuations are characterized 
by the Matsubara bosonic frequencies v = 27rnT. 

The resulting expression for the partition function may 
be simplified since the terms that survive the disorder- 
average in the sum of Eg. (|A12p are those for which^i 
qi = q2. Following Ref. [3l|, we disorder-average over the 
exponent in Eq. (|A6p . rather than over the free energy, 
to obtain an answer which is correct to leading order in 
(/iT+)~^. Finally we trace over the product of the 2x2 
matrices in Eq. (jA12p and integrate over A in Eq. (jA6p . 
In this way we obtain the partition function, Eq. ([7]). 
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